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Abstract. We have used measurements by the 
Fermi Large Area Telescope (LAT) to study the 
interstellar gamma-ray emission in a region of 
the second Galactic quadrant, at 100° < I < 145° 
and —15° < b < 30°. This region encompasses the 
prominent Gould-Belt clouds of Cassiopeia, Cepheus, 
and the Polaris flare, as well as large atomic and 
molecular complexes at larger distances in several 
spiral arms. The good kinematic separation in ve- 
locity between the local, Perseus, and outer arms, 
and the presence of massive complexes in each, make 
this region very well suited to probe the gamma-ray 
emission from the interstellar medium beyond the so- 
lar circle. The unprecedented sensitivity and angular 
resolution of the LAT provide improved constraints 
on the gradient of the cosmic-ray densities and on the 
increase of the CO-to-H 2 conversion factor, X C o, in 
the outer Galaxy. 
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I. Introduction 

Galactic interstellar gamma-ray emission is produced 
through the interactions of cosmic rays (CR) with the 
interstellar gas (via ir° production and Bremsstrahlung) 
and with the interstellar radiation field (via Inverse 
Compton interactions). Thus gamma rays can be used 
to trace cosmic rays in the Galaxy and the properties of 
the interstellar medium (ISM), in particular the total gas 
column densities. 

The interpretation of the observed gamma-ray diffuse 
emission is often based on two radio tracers of the in- 
terstellar gas. The 21 cm line of the hyperfine transition 
of atomic hydrogen Hi is used to derive its column 
density iV(Hi). The integrated 2.6 mm line intensity 
of the rotational transition of CO, Wco, is assumed 
to be proportional to the column density of molecular 
hydrogen N(H.2) through the conversion factor Xco, 
Xqo = N(R 2 )/Wco- The estimate of the Xco ratio 
is necessary for measuring H2 cloud masses from radio 
observations and it also impacts the derivation of the 
distribution of Galactic cosmic -ray sources from gamma- 
ray observations 11201 . but our current knowledge is 
limited, especially beyond the solar circle. Similarly, the 
gradient of the gamma-ray emissivity of the more diffuse 



Hi gas provides useful constraints on the gradient of 
the cosmic-ray densities in the outer Galaxy, but its 
derivation was so far limited by the sensitivity and 
angular resolution of gamma-ray telescopes. 

The situation has greatly improved on both grounds 
with the successful launch of the Fermi Gamma-ray 
Space Telescope on June 11, 2008. The Large Area 
Telescope (LAT) on board the Fermi mission [6| detects 
photons from 20 MeV to > 300 GeV. The sensitivity 
of the LAT is more than an order of magnitude greater 
with respect to the previous EGRET telescope on board 
the Compton Gamma-Ray Observatory; in addition, it 
has a superior single photon angular resolution (e.g., the 
68% containement angle for 1 GeV photons at normal 
incidence reaches 0.6° for the LAT with respect to 1.3° 
for EGRET). 

Here we present the analysis of the diffuse gamma-ray 
emission observed by the Fermi LAT in a selected region 
of the second Galactic quadrant, at 100° < I < 145° and 
— 15° < b < 30°, during the first 9 months of the science 
phase of the mission. 

II. Modelling of gamma-ray emission 

Our aim is to separate the gamma-ray emission from 
four regions along the line-of-sight in the second Galac- 
tic quadrant: 

1) the very nearby gas complexes in the Gould Belt, 
at hundreds of parsecs from the solar system [14J; 

2) the main part of the local arm, typically 1 kpc 

away; 

3) the Perseus arm, about 4 kpc away; 

4) the outer arm and beyond. 

In the (l,b) interval we have chosen the division is 
possible due to the good kinematic separation of the 
Hi and CO lines along the lines of sight and the low 
spatial degeneracy between the resulting maps in each 
region and gas phase. 

To estimate the gas column densities in each region 
we derived maps from the LAB survey [12| for Hi and 
from the moment-masked composite survey by Dame 
et al. [7 1 for CO. The following separation scheme 
was developed and applied to separate the four regions 
in velocity. The preparation of the gas maps started 
from preliminary velocity boundaries given in terms of 
Galactocentric annuli (R < 8.8 kpc, 8.8 kpc < R < 10 
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kpc, 10 kpc < R < 14 kpc, and R > 14 kpc) using 
a flat rotation curve with Rq = 8.5 kpc and a rotation 
velocity of 220 km s _1 . The velocity boundaries were 
adjusted for each line of sight to better separate the 
different cloud structures based on their coherence in 
the (l,b,v) phase space. These boundaries were used 
to produce iV(H i) and Wco maps in each region. The 
broad H i lines can spill over from one velocity band into 
the next. To correct for this cross-contamination between 
adjacent regions, each Hi spectrum was fitted by a 
combination of Gaussian lines. The latter were used to 
correct the integral in a particular velocity band from the 
spill-over from adjacent lines. The Hi column densities 
were derived applying the optical depth correction under 
the usual approximation of a uniform spin temperature 
T s = 125 K. 

An additional map was prepared to account for the so- 
called dark gas, i. e. neutral interstellar gas not traced 
by H i or CO. Grenier et al. iflOl reported a considerable 
amount of this gas at the interface between the phases 
traced by Hi and CO. The dark-gas map was derived 
from the E[(B - V)] map by Schlegel et al. [17], which 
traces the total dust column density in our Galaxy. This 
reddening map was fitted with a linear combination 
of the AT(Hi) and Wco maps previously described, 
to subtract the parts linearly correlated with this set 
of maps. The residual E[(B — V)] map traces the 
distribution of the gas column density not traced by H i 
and CO. 

The Inverse Compton (IC) emission expected over this 
narrow part of the sky is rather uniform and so indis- 
tinguishable from the isotropic background (the mixture 
of truly isotropic gamma-ray emission, unresolved point 
sources and residual background from misclassified in- 
teractions of cosmic rays in the LAT). Therefore we 
inserted only an isotropic term in the analysis model, 
to account for all of them. 

In this preliminary analysis we neglect any contribu- 
tion from the ionized hydrogen Hn. Thus, assuming 
that the ISM is transparent to gamma rays, that the 
characteristic diffusion lengths for cosmic electrons and 
protons exceed the dimensions of cloud complexes, and 
that cosmic rays penetrate clouds uniformly to their core, 
we can model the gamma-ray flux $ in a direction (I, b) 
by equation [T] 

$(l,b) = 5^[A-JV(Hi)(/,6),+B t -Wc3o(/,6).] + 
+ C + J2 D j-S (2) {l-lj,b-b J ) + 

+ E-E[(B-V)] ies (l,b) (1) 

The sum over i represents the combination of the four 
different Galactic regions. The parameters are the emis- 
sivity of Hi, A l7 and the emissivity per unit of CO 
emission intensity, B l . The isotropic term, C, accounts 
for all the components discussed above, including IC 
emission. The contribution from the 9 point sources 



found in the region in the LAT Bright Source List [2| 
at a significance level above 10 a is represented by the 
sum over j. For point sources we used the published 
positions, power-law spectra were fitted independently in 
each energy range. Three other point sources within 5° 
outside of the region boundaries were included, keeping 
their spectra as determined in the 6 month catalog 
internally available to the LAT Collaboration. The E 
parameter is the emissivity per unit of E[(B — V)] res . 

III. Gamma-ray analysis 

The gamma-ray data were obtained by the Fermi 
LAT during the period August 4 2008 - April 26 2009. 
The standard event selection was applied, selecting the 
Diffuse event class (the class with the least residual 
instrumental background), and rejecting events seen at a 
zenith angle greater than 100° and time intervals where 
the rocking angle respecting to zenith was greater than 
39°, in order to limit the contamination from Earth 
albedo gamma rays. 

The LAT data were analysed using a binned maximum 
Likelihood procedure with Poisson statistics, on a spatial 
grid with 0.5° spacing using the LAT Science Tool^. 
The analysis was done for four contiguous energy bands: 
0.3 - 0.6 GeV, 0.6 - 1 GeV, 1 - 2 GeV and 2-10 
GeV. A power-law spectrum was fitted independently for 
each diffuse component in each energy range, allowing a 
careful calculation of the exposures. Over such a narrow 
region of the sky we were able to perform a full con- 
volution with the energy-dependent PSE The analysis 
used the P6_V3 Instrument Response Functions; with 
respect to pre-flight Response [6| they include inflight 
calibration and account for the loss of efficiency due to 
pile-up events in the LAT ifTBI . 

We verified that the use of an IC emission map derived 
from the GALPROP cosmic -ray propagation code lfT9l . 
071 does not impact the resulting gas emissivities. The 
effect of the point sources was assessed by adding one 
by one the sources in order of decreasing detection 
significance: only the three brightest have an impact on 
the gas emissivities, except for the outer arm. Due to 
the small amount of gas in this region and the lack 
of significant features in its spatial distribution, the gas 
emissivity responds to the presence of point sources 
along the Galactic plane. Thus the current A4 value must 
be regarded only as an upper limit and will not be used 
in the following discussion. 

A detailed description of the analysis, including a 
careful treatment of the systematics, will be given in 
a paper in preparation within the Fermi LAT Col- 
laboration. In the following section we discuss some 
preliminary results. 

IV. Results 

A. H 1 emissivities 

The integrated Hi emissivities measured in the four 
energy bands are given in Fig. [JJ The present values 

'http://fermi.gsfc.nasa.gov/ssc/data/analysis/documentation/Cicerone/ 
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Fig. 1. (PRELIMINARY) Dependence on Galactocentric radius R of the Hi emissivity measured in the four energy bands: a) 0.3 - 0.6 GeV, 
b) 0.6 - 1 GeV, c) 1 - 2 GeV, d) 2 - 10 GeV. Black dashed lines give the predictions from the GALPROP model 54_77Xvarh7S, i.e. the 
sum of the gas emissivity due to 7r° decay and Bremsstrahlung, scaled by +15%. The three regions (blue crosses) spanning different ranges in 
R encompass the Gould Belt, the local arm and the Perseus arm (from left to right); the error bars include systematics due to the uncertainties 
on the event selection efficiency. 



in the Gould Belt are consistent within 10% with the 
local H i emissivity spectrum derived from LAT data in 
a mid-latitude region of the third Galactic quadrant |4j. 

The physical interpretation of the emissivities is based 
on the comparison with a GALPROP model. GALPROP 
is a numerical code for cosmic -ray propagation in the 
Galaxy, see e.g. fl9l . ET1 . widely used to analyse 
LAT data. The program has recently been updated 
with the inclusion of new astronomical surveys, more 
accurate description of the interstellar gas, improvements 
in the modelling of physical processes, the cosmic- 
ray electron spectrum derived from Fermi data 0] and 
tuning to the diffuse gamma-ray emission seen by Fermi 
itself E2. The GALPROP run used for this work, 
5 4_7 7Xvarh7S, is consistent with directly measured 
CR spectra and the non-confirmation of the EGRET GeV 
excess 0, lfl5l . It has been scaled by +15% in Fig. Q]to 
provide a straightforward comparison of the emissivity 
decrease towards the outer Galaxy. This factor is the 
average of the data-over-GALPROP ratios found in the 
four energy bands for the Gould Belt. 

Provided this small intensity normalization, the con- 
sistency between the data and GALPROP points in the 
Gould Belt at all energies shows that the gamma-ray 



emissivity of local H i can be explained by the cosmic- 
ray spectrum measured at Earth. This confirms for very 
nearby gas the conclusion reached from the analysis of 
LAT data over broader regions of the sky Q, fl5l : LAT 
observations are inconsistent with the GeV excess seen 
by EGRET over the whole sky and in specific local H i 
complexes |9|. No significant variations in emissivity are 
found between the Gould Belt and the local arm, as can 
be expected from the large coupling length of ~ 1.75 
kpc found between cosmic protons and interstellar gas 

mi. 

The emissivity profiles in Fig. Q] significantly decrease 
with Galactocentric distance, as expected from the dis- 
tribution of candidate cosmic-ray sources used in the 
GALPROP model. However, the measured decrease is 
found to be systematically smaller from the local to the 
Perseus arm than the GALPROP prediction, especially at 
high energies where the higher angular resolution allows 
a better component separation (Fig. [Tc1and[Tdb. 

B. Calibration of molecular masses 

The emissivities can be used to estimate the CO-to- 
H2 conversion factor. In each energy range and each 
region we can derive Xqo% — B l /2A l . Derived as a flux 
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ratio between two components with the same spectral 
shape, the Xco factor is not affected by systematic 
uncertainties in the instrument efficiency. It is also 
insensitive to the presence of isolated point sources not 
included in the analysis, because the large angular extent 
and characteristic shapes of the clouds provide better 
discrimination in the fitting procedure, thanks to the 
good angular resolution of the LAT. However, the Xco 
factor is sensitive to the inclusion of the dark gas in 
the model. The impact is not highly significant in the 
very nearby Gould-Belt clouds, where we have ample 
resolution and little degeneracy between the Hi, dark, 
and CO phases. But, since the three phases blend more 
when seen from a larger distance, the separation is less 
effective and the Xco measurement in the more distant 
complexes becomes more sensitive to the inclusion of 
the dark gas. For instance, the Xco estimate in the local 
and Perseus arms can increase by 30% if the dark gas 
is not included in the analysis. The estimate of Xco 
is also sensitive to the level of degeneracy between 
N(H i) and Wco maps that is inherent to the physical 
structure of the clouds. Since the angular resolution of 
the LAT is strongly energy-dependent, this can lead to an 
energy-dependence of the Xco measurement for distant 
clouds. The small spread of Xco measurements with 
energy can be considered as a systematic uncertainty on 
the component separation. The results are summarized 
in Fig. 12 where Xco is evaluated for each region as 
the weighted average of the values obtained in the four 
energy ranges and the error bars include these systematic 
uncertainties. 

The value measured in the Gould-Belt clouds agrees 
with the results obtained by EGRET for the Cepheus and 
Polaris flares J8), and is consistent with estimates for 
other clouds in the solar neighborhood [9|. The Gould- 
Belt clouds are within 300 pc from the solar system, so, 
given the < 0.6° angular resolution of the LAT above 
1 GeV, they are mapped in gamma-rays with a spatial 
accuracy < 3 pc. They also show a very low attitude 
at forming massive stars. Therefore it is very unlikely 
that the corresponding Xco measurement is biased by a 
population of unresolved sources closely following the 
structure of the clouds. This bias is unlikely too in the 
local-arm region, but source contamination cannot be 
ruled out for the Xco measurement in the more remote 
Perseus arm. 

The present results suggest an increase of Xco m 
the outer Galaxy, as already deduced by many other 
observations (see e.g. Sodroski et al. |[T8l and the other 
references cited in ll20"l ). Whether the present Xco gra- 
dient can be fully attributed to the metallicity gradient, 
or partially to unresolved sources, H i and CO separation 
problems, or gas not traced by H i and CO, needs further 
investigation. In any case, the results shown in Fig. [2] 
indicate smaller values of Xco with respect to those 
used by Strong et al. |20) in GALPROP: the present 
measurement is significantly lower in the outer Galaxy. 
The values preliminarily obtained from LAT data are 




Fig. 2. (PRELIMINARY) Variation in function of Galactocentric ra- 
dius of the Xqq conversion factor. Blue crosses are our measurements 
in the distance ranges corresponding to the Gould Belt, the local arm 
and the Perseus arm (from left to right). Dashed black lines represent 
the values used by Strong et al. [20] in GALPROP. The solid red line is 
the conversion function determined for the outer Galaxy by Nakanishi 
& Sofue [131 from CO data and virial masses (adapted to the rotation 
curve assumed for our analysis). 



also systematically smaller than the Xco(R) relation 
determined from CO data and virial masses beyond the 
solar circle ll5l. lfT3l. 
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